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Abstract 

The phase transitions occurring in the frustrated Ising square antiferromagnet 
with Erst- (Ji < 0) and second- (J 2 < 0) nearest-neighbor interactions are studied 
within the framework of the effective-held theory with correlations based on the 
different cluster sizes and for a wide range of i? = J 2 /1 Ji |. Despite the simplicity of 
the model, it has proved difficult to precisely determine the order of the phase tran¬ 
sitions. In contrast to the previous effective-held study, we have found a hrst-order 
transition line in the region close to R = —0.5 not only between the superantifer¬ 
romagnetic and paramagnetic {R < —0.5) but also between antiferromagnetic and 
paramagnetic {R > —0.5) phases. 
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1. Introduction 


The Ising square lattice with nearest-neighbor (nn) interactions (Ji) is a rare 
instance of the exactly soluble model [1, 2] which shows a phase transition. However, 
adding next-near-neighbor (nnn) interactions (J 2 ) or a magnetic held (or both), 
makes the model no longer exactly soluble and only some approximate solutions 
are possible to attack this more general problem. The behavior of this model is 
of interest, not only because of the existing theoretical questions, but also because 
the quasi-two-dimensional anisotropic antiferromagnet like K 2 C 0 F 4 compound [3] 
can be represented by Ising model at least within a fair approximation. Since the 
existence of nnn interactions and the magnetic held may give rise to other phases 
with diherent types of phase transitions and multicritical points [4], we restrict 
our discussions of this model to the antiferromagnetic system with zero magnetic 
held. In this case a nnn antiferromagnetic interaction represents the simplest way 
to incorporate frustration in the square bipartite lattice. The Hamiltonian of the 
model with competing antiferromagnetic interactions between nn (Ji < 0 ) and nnn 
(J 2 < 0 ) spins is given by 

H J\ ^ ^ SiSj J 2 ^ ^ (1) 

(*J> (*>*2> 

where s* = ±1 and the hrst summation is carried out only over nn pairs of spins 
and the sums in the second term run over nnn pairs of spins. We note that in the 
zero magnetic held under appropriate transformations, the antiferromagnetic and 
the ferromagnetic states are equivalent. Therefore, the results are independ of the 
sign of Ji and may be compared with those obtained earlier for the ferromagnetic 
case (Ji > 0 ). 

The ground state of ([T]) is well known. Namely, for the value of the frustration pa¬ 
rameter R = J 2 /\Ji\ > —0.5 the ground state is the Neel antiferromagnet (AF) and 
in the case oi R < —0.5 the system orders in alternate ferromagnetic rows (columns) 
of opposite oriented spins (superantiferromagnet) [5]. The critical point separating 
these two phases lies at i? = —0.5, where the transition temperature is suppressed 
to zero. At this point the ground-state is highly degenerate, a phenomenon which 
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is responsible for the vanishing of long-range order at hnite temperatures, and the 
system behaves like a one-dimensional Ising model [6 — 9]. 

Even though the investigations of the frustrated Ji — J 2 Ising model dehned by 
the Hamiltonian ([1]) have a long history, its hnite-temperature phase transition re¬ 
mains controversial until now. Indeed, the early numerical and analytic approaches 
supported the idea that the transition is always continuous for R < —0.5, but 
with critical exponents that vary with R (weak universality) [4, 6,10,11]. However, 
some approximate studies [8, 9,12] and recent Monte Carlo method [13,14] have 
found a hrst-order transition for —1.1 ^ R < —0.5. On the other hand, very re¬ 
cent cluster mean-held calculation [15] with cluster of size 4x4 also gives change 
the antiferromagnetic-paramagnetic transition to hrst-order in a narrow region of 
R > —0.5. Unfortunately, there are no clear indications from previous Monte Carlo 
simulations of the transition being hrst-order in this regime. Therefore, it may be 
worthwhile to examine this issue more systematically by using diherent size of the 
clusters in order to carefully analyze the ehect of the cluster size on the occurrence of 
hrst-order transitions. To this end we study the Ji — J 2 model in its parameter space 
using an ehective-held theory based on two-, four-, six-, and nine-spin clusters. This 
approach is based on the diherential operator technique introduced into exact Ising 
spin identities and has been successfully applied to a variety of spin-| and higher 
spin problems (for a review see, e.g.. Ref. [16]) including a geometrically frustrated 
triangular lattice Ising AF [17 — 19]. It is important that the present ehective-held 
theory allows us to treat large clusters in a simpler and more efficient computational 
manner. Of course, the ehective-held calculation cannot give any information on the 
true critical exponents, although the second-neighbor interactions should in general 
lead to critical exponents which diher from those of the nearest-neighbor square 
Ising lattice [20]. 

2. Theory 

The starting point for the statistics of an Ising spin system, which is much more 
useful for our purpose, is the following generalized Callen-Suzuki [21, 22] exact iden- 
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tity 


{0{n}) = 


Tr{n} [ 0{n} exp {-(3H{n})] 


( 2 ) 


Tr{n}[exp(-/3i/{„})] 

where the partial trace Tr is to be taken over the set {n} of spin variables specified 
by the clnster spin Hamiltonian H^n}- Here 0{„} denotes any arbitrary spin fnnction 
inclnding the set of all {n} spin variables (finite cluster) and (• • •) denotes the usual 
thermal average. 

In order to take into account effects of frustration within the present effective- 
field theory, it is necessary to consider at least a two-spin cluster. In this approach, 
we select two nn spins, labeled i and j, which interact with other nn and nnn spins 
from the neighborhood [23]. Hence, the multi-spin Hamiltonian for the AF 
two-spin cluster (n = 2) on the square lattice (Fig. 1(a)) is given by 



- - sfhf, 

(3) 



(4) 


with 


21 = 1 22 = 1 J'l = l ^2 = 1 

where sf is a spin variable on sublattice A and is a spin variable on sublattice B, 
the superscript AF denotes the antiferromagnetic system, and the terms A = j and 
ji = i are excluded from summations over A and ji, respectively. Now, according 
the exact Callen-Suzuki identity ([2]), the normalized staggered magnetization ttiaf 
associated with the AF two-spin cluster is given by 


rriAF = 


A 


Sf) = 


D. 


" )fAF{x,y) 


x=0,y=0 


(5) 


where = d/dy (/i = x,y) are the differential operators and function fAF{.x,y) is 
defined by 

f ( ^ sinhA(a:-|/) 

fAF{x,y) = —( 6 ) 


cosh A (a: — y) + cosh/3 (a: -|- y) 

At this point one should notice that the neighborhood of the sites i and j of the 
two-spin cluster for the Ji — J 2 model on a square lattice contain a set of common 
spins, namely the spins at the sites labeled by (A, 22 ) or (ji. A) in Fig- 1(a)- These 
spins are frustrated directly within two-spin cluster theory, which is not the case 
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of the one-spin cluster approximation. Now, taking this into account and assuming 
the statistical independence of lattice sites, Eq. (|5]) may be rewritten as 


rriAF = 


^a:(l)^y(2) -|- Bx(l)By(2) + UIb ( Ax{l)By(2) -|- Ay(2)Bx{l) 


X 


+ By{l)Bx{2) -|- niA [ + Ax(2)By{l) 


X A^{1) + mBB^{l) Ay{l) + mABy{l) 


X 


A^{2) + mABa,{2) ) ( Ay{2) mBByi2) ) fAF{x, y) 


a;=0,y=0 


(7) 


where A^(z/) = cosh( B^iy) = sinh( (z/ = 1,2), and = (s“) 
(a = A, B) are the sublattice magnetizations per site. This equation of state is quite 
superior to that obtained from the standard mean-held theory, since in the present 
framework relations like {Sg)‘^^ = 1 and = Sg, for all r, are taken exactly 

into account through the van der Waerden identity, exp(as“) = cosh(a) -|-s"sinh(a), 
while in the usual molecular-held theory all the self- and multi-spin correlations 
are neglected. Since for the AF A — J 2 model at zero magnetic held we have 
mAF = 1 x 1 A = Eq. ([7]) can be hnally recast in the form 

4 

rriAF = ^m^af^ , ( 8 ) 

n=0 

where the coefficients which depend on T and R, can easily be calculated 

within the symbolic programming by using the mathematical relation exp(AZl 3 ; -|- 
lDy)fAF{,x, y) = fAF{,x+\, 2 /-I- 7 ). Because the hnal expressions for these coefficients 
are lengthy, their explicit form is omitted. We also note that in obtaining Eq. (IH]) we 
have made use of the fact that fAF{x,y) = —fAF{—x, —y) and therefore only odd 
diherential operator functions give nonzero contributions. 

In order to determine the phase diagram of the AF A — J 2 model, we should solve 
Eq. ([H]) for a given value of the frustration parameter R and look for the temperature 
at which the magnetization (order parameter) itlaf goes to zero. However, for some 
values of i?, the order parameter goes to zero discontinuously, i. e., the transition 
becomes hrst order. To analyze hrst-order transitions, one needs to calculate the free 
energy for the AF and paramagnetic (P) phases and to hnd a point of intersection. 
Because the expression for the free energy in this effective-held theory does not exist. 
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it will be extrapolated with the help of the relation for the equilibrium value of the 
order parameter ([8]) as follows [24]: 



(9) 


n=0 


where Fq{T,R) is the free energy of the disordered (paramagnetic) phase and m is 


the order parameter which takes the value niAF at thermodynamic equilibrium. We 


note that relation (|9]) corresponds to a Landau free energy expansion in the order 
parameter truncated in the term. 

Thus the equilibrium staggered magnetization is the value of the order parameter 
which minimizes the free energy given by Eq. (|9]). Using the equilibrium condition 


dF^^{T,R,m) 


( 10 ) 


0 . 


dm 


we recover Eq. (jHj) for the equilibrium staggered magnetization. Then a critical 
temperature and a tricritical point, at which the phase transition changes from 
second order to hrst order, are determined by the following conditions [25]: (i) the 
second-order transition line when l — K^^ = 0 and < 0, and (ii) the tricritical 

point (TCP) when 1 — = 0, = 0, if < 0. However, the hrst- 

order phase transition line is evaluated by solving simultaneously two transcendental 
equations, namely the equilibrium condition fITII]) and the equation F^^ {T, R,m) = 
Ff){T, R), which corresponds to the point of intersection of the free energies for the 
AF and P phases. It is worth noticing that if we use the above given relations to 
obtain the critical and tricritical points, the results will coincide with those obtained 
from Eq. ([8|). As has been pointed out in [25], these results justify our procedure. We 
also note that a similar methodology of obtaining the free energy of the model within 
the effective-held theory, employing the two-spin cluster, has been proposed in Ref. 
[26] had used recently to investigate phase diagram of the frustrated Heisenberg 
model [27]. 

On the other hand, the choice of a two-spin cluster for the superantiferromagnet 
(SAF) is not unambiguous. Indeed, one can chose a cluster with two spins in the 
vertical or horizontal direction. As we will see below, only results obtained for 
the latter case are consistent with the ground-state behavior of the model. The 
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Hamiltonian for such a choice of the cluster with two spins in the horizontal direction 
(see the scheme in Fig. 1(b)) is given by 


Tj SAF — _ j 0 ^ 0 ^ — Q^h‘ 
'tiSj Sj Sj ti. 


A„A ^Ai^SAF ^AuSAF 
i 


s^h^- 


( 11 ) 


where the corresponding ‘helds’ and are given by 


^ sf) + ^2 ^ sf^, - Ji(sf^=i + gu )+ -^2 


2l =2 


22 = 1 


ii=2 




i 2 =i 


( 12 ) 

Then, in the same way as for the AF, one can derive equation of state for the SAF 
in the form 


msAF = 


Ax(l)Ay(2) + Bx{l)By(2) + rriB (A3;(l)i?y(2) + Ay(2)Bx{l) 


X 


Aj,(l)Aa;(2) + By{l)Bx{2) + ruB [ Ay{l)Bx{2) + A^(2)By{l) 


X ^ 2 ,( 1 ) + mAB^{l) Ay{l) + mABy{l) 


X 


A^{2) + msHa,(2)) ( Ay(2) + mBBy{2)) fsAF^x, y) 


x= 0 ,y =0 


.( 13 ) 


where ttisaf = (l(s(^ + sf)) = ruA = —'m-s and function fsAF^x, y) is given by 


fsAF{x,y) = 


sinh/S(a; + y) 


cosh/9(a; + y) + cosh /3{x — y) 


(14) 


Further, using the same procedure as above for the AF, we derive the free energy 
functional to determine numerically the transition lines, including the TCP between 
the SAF and P phases. 

One can continue this series of approximations to consider larger and larger 
clusters and as a consequence, better results are expected. Therefore, the next 
step for the square lattice is a four-spin cluster (n = 4) containing the symme¬ 
try of both the AF and SAF phases. In this case, we start from Eq. ([2]) with 
0{n} = {sf + sf — sf — sf )/4 for the AF and C{„} = (sf -1- sf — sf — sf )/4 for the 
SAF. The formalism is developed in the same way as above for the two-spin cluster. 
For instance, one derives for the AF the equation analogous to ([5]), which reads 


rriAF 


l^^hf^Di+h^‘^D2+hf^D3+hff'D4 


fAF{Xl,X2,X3,Xi) 


{xi=0} 


(15) 
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where each variable Xi{i = 1,2, 3,4) of the function fAF{xi,X 2 ,X 3 ,X 4 ) (which, for 
brevity, is not presented explicitly here) is associated with a differential operator 
Di = d/dxi{i = 1,2, 3,4). The corresponding ‘helds’ and hf^ (see 

also [12]) are given by 


k 


AF 


h 


AF 


2 3 2 3 


"^1 4 

+ 


ii=l 

^2 = 1 

ii=i i2=i 

2 

3 

2 3 


+ J2^ 

4^ = 4 + si 

ki=l 

k2=l 

/l = l /2 = 1 


(16) 


Then, using of the same procedure as for the two-spin cluster, one determines the 
second- and first-order transition lines inclnding the TCP between the AF and P 
phases. Finally, to get a more convincing evidence of the existence of a TCP in the 
phase diagram, we also consider six- and nine-spin clusters. In this case, the clusters 
consist of two or fonr sqnare plaqnettes, respectively, and contain more information 
abont the lattice topology. However, the analytical calculations for such large clus¬ 
ters would have been very lengthy and tedious, therefore, the results were obtained 
in a completely nnmerical way within the symbolic programming by using Mathe- 
matica software package. It should be noted here that the calculation times for the 
large clusters become rather long even using the symbolic programming. Therefore, 
the highest approximation used to study the frustrated Ising square lattice is the 
one based on the nine-spin cluster. 


3. Results 


The numerical results for the critical temperature /cbTat/I Ji| versus R for various 
cluster sizes are shown in Figs. 2, 3, and 5. In these hgnres the solid lines indicate 
the second-order phase transitions, while the dashed lines represent the hrst-order 
ones. The black circles denote the positions of TCPs at which the phase transitions 
change from second to first order ones. 

First, from all the hgnres it can be seen that the hrst-order transition temper- 
atnres between the AF and P phases as fnnctions of the frustration parameter 
R, for the considered cluster sizes, approaches zero when R = —0.5, in agreement 



with the ground-state result. For the nonfrustrated model (i? = 0), we hud that 
values of /cbTat/IJiI are 3.0250,2.9197,2.8759, and 2.8185 for the two-, four-, six-, 
and nine-spin clusters, respectively, which shows a relatively slow convergence to the 
exact value of A;bT/v/|Ji| = 2.2692 with increasing cluster size. Further, it is seen 
that the critical temperature decreases gradually to the TCP, when the frustration 
parameter approaches some negative value Rt- Our estimates for the coordinates of 
the TCP (M/I Jil; are (1.3720; -0.2973), (1.3066; -0.3070), (1.2329; -0.3240), 
and (1.1894; —0.3340) for the two-, four-, six-, and nine-spin cluster approximations, 
respectively. Thus, within the present approach the tricritical temperature decreases 
with increasing cluster size. However, the relative decrease of the tricritical temper¬ 
ature becomes gradually smaller with the cluster size. For instance, the tricritical 
temperatnre for the four-spin cluster is about 5.00% lower than that for the two- 
spin cluster but for the nine-spin cluster the derease is only 3.66% in comparison 
with the six-spin cluster. These results indicate that there could be a narrow region 
of hrst-order transitions for R > —0.5 and that the cluster-size has a considerable 
effect on the existence and location of the TCP at which the phase transition be¬ 
tween the AF and P phases changes from second order to first one. It should be 
noted here that the tricritical behavior in this region has also been predicted by 
a cluster mean-held approach based on the 4x4 cluster but not for the 2x2 cluster 
[15]. Such behavior of the model for R > —0.5 is also supported by the fact that 
at i? = —0.5 there is, as discussed above, a hrst-order transition at T = 0 between 
the AF and SAF phases. We believe that the hrst-order behavior of the model in 
the region R > —0.5 cannot be completely ruled out and should be examined close 
to i? = —0.5 not only by various approximate theories but also by more reliable 
techniques, such as e.g. Monte Carlo approach, to dehnitely conhrm or to rule out 
this scenario. Note that this model has also been investigated by the ehective-held 
theory based on the four-spin cluster in [12]. Unfortnnately, some of the nnmerical 
resnlts obtained in that work are defective, and the anthors overlooked the existence 
of the TCP on the transition line between the ferromagnetic (or antiferromagnetic) 
and paramagnetic phases. It is worth mentioning that we also checked numerically 
that the obtained phase diagrams are independent of the sign of A. Of course, for 
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Ji > 0 (ferromagnetic case), the order parameter mp must be defined as the direct 
(not staggered) magnetization. 

On the other hand, it is generally accepted now that the transition line between 
the SAF and P phases is a second order for the full range of R < —1.1 and a hrst 
order for —1.1 ^ R < —0.5. Such behavior is conhrmed not only by approximate 
methods [8, 9,12], but also by recent Monte Carlo studies [13 — 15, 28]. We note that 
the region of —1.1 < i? < —0.67 turns out to be very difficult to study by numerical 
means, nevertheless, the results of the recent study [28] indicate only pseudo-first 
order behavior in this region but the true nature should be second order. Our 
effective-field results obtained for various cluster sizes support this picture. How¬ 
ever, as seen from Fig. 2, the transition line as well as the coordinates of the TCP 
obtained within the effective-held theory based on the two-spin cluster for R < —0.5 
depend on the selection of the two neighboring spins. Namely, if we chose a cluster 
with two spins in the horizontal direction (see the scheme in Fig. 1(b)), the coordi¬ 
nates of the TCP are (2.7178; —1.1098), while if consider a cluster with two spins in 
the vertical direction (the scheme similar to that of Fig. 1(b) is not shown here), the 
coordinates of the TCP are (2.2594; —0.9991). Moreover, the hrst-order transition 
for the latter case terminates at i? = —0.6168 (see the inset of Fig. 2), which is 
inconsistent with the ground-state behavior of the model. Therefore, the effective- 
held theory based on the cluster with two spins chosen in the vertical direction is 
not appropriate for the SAF state due to its incompatibility with the symmetry of 
the SAF state. 

It is evident, that this lack of symmetry for the SAF state is absent in the four- 
or nine-spin clusters. Consequently, one can hnd that the corresponding hrst-order 
phase transitions terminate at i? = —0.5, as expected from the ground-state ar¬ 
guments. Since the phase diagrams for these clusters are qualitatively the same, 
in Fig. 3 we show only the phase diagram for the nine-spin cluster. Note also 
that the TCPs for the four- and nine-spin clusters are located uniqly and their 
coordinates are (2.4620;—1.0387) and (2.2565;—0.9748), respectively. The values 
Rt = —1.0387 and Rt = —0.9748 may be compared to those of the cluster-variation 
method {Rt ~ —1.144) [8,9], Monte Carlo study {Rt ~ —0.9) [14], and recent 
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Monte Carlo approach [Rt ~ —0.67) [15,28]. On the other hand, for the effective- 
field theory based on the six-spin cluster there are two possibilities how to select 
the cluster for the SAF, namely in the horizontal or in the vertical direction (see 
Fig. 4). It is clear now that the clusters consist of different number of spins on sub¬ 
lattices A and B. Consequently, the coordinates of the TCP are (2.4018; —1.0111) 
or (2.2194; —0.9666), respectively, and the corresponding phase diagrams are shown 
in Fig. 5. But now both clusters consist of two square plaquettes, which are sym¬ 
metrical with respect to the SAF state. Therefore, the phase diagrams for these 
clusters differ only quantitatively and end at R = —0.5, which is exact. 

In order to conhrm the prediction of the first- and second-order phase transitions, 
let us examine temperature dependencies of the order parameters uiaf and msAF 
for the nine-spin cluster, when the value of R is changed. As seen from Fig. 6(a), the 
order parameter niAF falls smoothly to zero when temperature increases from zero 
to ksTAf/lJil, characterizing a second-order phase transition. Similarly, the ttiaf 
also reduced to zero continuously at the TCP (see curve labeled —0.3340). On the 
other hand, below the TCP, the stable solution of tuaf becomes discontinuous at the 
first-order phase transition and this discontinuity increases with R going to —0.5. 
The curves for R = —0.35 and R = —0.4 are examples of such behavior, where the 
hrst-order transition is indicated by a vertical dashed line. A qualitatively similar 
behavior exhibits the order parameter fusaf (Fig. 6(b)). 

4. Conclusions 

In this work, we have investigated the effects of the frustration parameter R = 
J 2 /\Ji\ on the phase diagram of the Ising antiferromagnet on the square lattice by 
the use of the effective-held theory with correlations and approximations based on 
diherent cluster sizes. For diherent hnite-clusters approximations it is found that 
the phase transitions between the AF and P phases {R > —0.5) close to R = —0.5 
are of hrst order. This behavior was conhrmed by solving the Hamiltonian using 
two-, four-, six-, and nine-spin clusters. Our results indicate that the cluster size 
plays an important role for determining the location of the TCP in the frustrated 
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square lattice. This behavior is consistent with that obtained within a variational 
mean-held theory [15], which predicts the TCP for the 4x4 cluster but not for the 
2 x2 cluster. Of course, these are the effective-held results, therefore, further Monte 
Carlo simulations or more reliable calculations in this region are called for. 

On the other hand, we have shown that the frustrated square lattice exhibits the 
TCP, at which the phase transition changes from second order to hrst order, on the 
line between SAP and P phases, in agreement with previous approximate as well as 
Monte Carlo studies. However, within present ehective-held theory, the position of 
the TCP for some clusters (two- and six-spin clusters) depends on the symmetry of 
the selected cluster. This behavior is due to a lower symmetry of the SAP phase in 
comparison to the AP phase. 
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Figure captions 


Figure 1: Ground-state configurations of the Ji — J 2 Ising model on the square 
lattice showing (a) aniferromagnetic and (b) superantiferromagnetic states for the 
two-spin cluster approximation dehned by the Hamiltonians ([3]) and flTT]) . respec¬ 
tively. Two sublattices are marked by black and white circles. 

Figure 2: Phase diagram in the coupling-temperature plane for the Ji — J 2 Ising 
model on the square lattice based on the two-spin cluster. Solid and dashed lines 
indicate second- and hrst-order transitions, respectively, while the black circles de¬ 
note the position of a tricritical point. AF and SAF are the antiferromagnetic 
and superantiferromagnetic ordered phases, respectively, and P is the paramagnetic 
phase. The inset shows phase diagram when the two-spin cluster is chosen in the 
vertical direction. 

Figure 3: Phase diagram in the coupling-temperature plane for the Ji — J 2 Ising 
model on the square lattice based on the nine-spin cluster. Solid and dashed lines in¬ 
dicate second- and hrst-order transitions, respectively, while the black circles denote 
the position of a tricritical point. AF and SAF are the antiferromagnetic and su¬ 
per antiferromagnetic ordered phases, respectively, and P is the paramagnetic phase. 

Figure 4: Two options of the six-spin cluster with spins Sj, s^, s^, s^, Sm, and for 
the superantiferromagnetic arrangement on the square lattice: (a) in the horizontal 
direction and (b) in the vertical direction. Two sublattices are marked by black and 
white circles. Each cluster has 14 nn and nnn spins. 

Figure 5: The same as in Fig. 2, but for the six-spin cluster. 

Figure 6: Temperature dependences of the order parameters of mAp in (a) and 
msAF in (b) for the Ji — J 2 Ising model on the square lattice based on the nine-spin 
cluster, when the frustration parameter R is changed. The dashed lines indicate the 
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first-order transitions. 
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